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ABSTRACT 


The  CH-46  Tie  Bar  is  a  multiple  leaf,  stainless  steel  system  that  attaches  the  rotor  blade  to 
the  rotating  hub  and  carries  rotor  blade  centrifugal  force.  The  Tie  Bar  twists  as  the  rotor  blade 
feathers  and  bends  slightly  as  the  Pitch  Housing  bends.  Originally  designed  to  last  3,000  hours 
without  any  component  breakages,  the  Tie  Bar  has  experienced  field  failures  that  have  necessitated 
inspections  every  10  flight  hours.  Traditional  fatigue  qualification  and  life  methodology  have  not 
reconciled  the  unusually  large  number  of  CH-46  Tie  Bar  component  failures.  One  CH-46 
experienced  15,  10,  and  6  component  failures  on  the  three  aft  rotor  head  Tie  Bars  after  only  93 
flight  hours.  One  complete  Tie  Bar  (140  components)  failed  and  caused  an  aircraft  crash. 

This  research  developed  a  probability  modeling  predictive  method  that  can  aid  the  design 
of  a  Tie  Bar  that  will  meet  desired  life  and  reliability  levels.  Laboratory  coupon  tests  were 
interpreted  to  provide  material  properties  for  several  manufacturing  processes.  Finite  element 
analysis  of  the  current  Tie  Bar  design  as  well  as  two  modified  designs  was  used  to  determine  the 
Tie  Bar  stress  state  for  many  component  failure  combinations.  The  reduced  constraint  design 
produced  a  dramatic  stress  concentration  reduction  and  may  provide  large  gains  in  life  over  the 
current  and  elliptical  slot  designs. 
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1.  INTRODUCTION 


A.  CH-46  HELICOPTER 

The  CH-46  helicopter,  pictured  below  in  Figure  I-l  and  Figure  1-2,  is  operated  by  both 
the  U.S.  Navy  and  the  U  S.  Marine  Corps.  Each  rotary  wing  head  on  the  CH-46  has  three  rotor 
blades.  A  diagram  of  the  rotary  wing  head  is  shown  in  Figure  1-3.  This  investigation  addresses  the 
analysis  and  prediction  of  field  failures  of  the  Tie  Bars  and  options  for  Tie  Bar  redesign. 


Figure  I-l  CH-46  Helicopter,  Side  View 
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B.  PITCH  ASSEMBLY 

The  exploded  view  of  the  CH-46  Rotary  Wing  Head  pitch  assembly.  Figure  1-4,  depicts 
the  Tie  Bar  and  the  other  rotor  blade  pitch  change  components.  The  cutaway  view.  Figure  1-5, 
illustrates  the  Tie  Bar’s  location  inside  the  Pitch  Housing  and  Pitch  Shaft. 

The  Tie  Bar  carries  the  rotor  blade  centrifugal  force.  One  end  of  the  Tie  Bar  is  connected 
to  the  Pitch  Housing;  the  other  end  is  connected  to  the  Pitch  Shaft.  There  are  roller  bearings 
between  the  pitch  shaft  and  pitch  housing  that  allow  relative  rotation  and  sliding  between  the  two 
parts.  Since  the  Tie  Bar  is  connected  to  both  the  Pitch  Housing  and  the  Pitch  Shaft  it  twists  as 
control  forces  are  applied  to  the  Pitch  Housing  to  feather  the  rotor  blade.  To  reduce  control  loads, 
the  Tie  Bar  is  designed  to  be  soft  in  torsion.  A  flatwise  bending  moment  due  to  flapping  applies  a 
bending  moment  to  the  pitch  shaft  and  pitch  housing  which  in  turn  bends  the  Tie  Bar. 
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Figure  1-4  CH-46  Rotary  Wing  Head  Pitch  Change  Assembly,  Exploded  View 
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Figure  1-5  CH-46  Rotary  Wing  Head  Pitch  Change  Assembly,  Cut-Away  View 
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C.  TIE  BAR 


The  Tie  Bar  is  constructed  of  301  Stainless  Steel.  It  consists  of  35  slotted  straps  that  are 
stacked  with  301  Stainless  Steel  spacers  in  between.  The  slots,  in  essence,  make  each  individual 
Tie  Bar  strap  a  redundant  four  component  assembly.  The  “components”  are  lettered  A  through  D 
in  Figure  1-6.  The  Tie  Bar  stack  is  held  together  by  six  hollow  pins  that  are  expanded  by  driving  a 
metal  ball  through  them  during  assembly. 


Figure  1-6  Tie  Bar  Strap  Geometry 

Designed  to  last  3,000  hours  without  any  component  breakages,  the  Tie  Bar  is  now 
inspected  every  10  flight  hours.  Normal  fatigue  qualification  and  life  methodology  have  been 
unsuccessful  in  predicting  an  unusually  large  number  of  CH-46  Tie  Bar  component  failures.  One 
CH-46  experienced  15, 10,  and  6  component  failures  on  the  three  aft  rotor  head  Tie  Bars  after  only 
93  flight  hours.  A  complete  Tie  Bar  failure  has  caused  one  aircraft  crash.  The  cause  for  early  in- 
service  failure  of  Tie  Bar  elements  has  not  been  determined.  Improved  manufacturing  methods  and 
quality  control  have  not  eliminated  the  problem. 

D.  PROBABILITY  MODELING 

To  predict  the  life  of  a  part  by  statistics  with  a  high  confidence  level  for  safety  related 
components,  the  number  of  test  specimens  must  be  at  least  an  order  of  magnitude  higher  than  the 
desired  reliability.  For  example,  to  obtain  a  reliability  level  of  one  failed  part  in  100,000  (10^),  at 
least  1  million  (10^^')  parts  need  to  be  tested.  Clearly,  this  is  impractical.  In  addition,  the  statistical 
approach  requires  building  the  actual  part  before  the  reliability  level  can  be  determined.  To  solve 
the  Tie  Bar  failure  problem,  a  modeling  method  needs  to  be  developed  that  can  aid  the  design  of  a 
Tie  Bar  that  will  meet  desired  life  and  reliability  levels. 


5 


Since  a  statistical  approach  is  unrealistic,  a  probability  approach  has  been  chosen. 
Through  a  combination  of  finite  element  modeling,  coupon  test  analysis,  fatigue  modeling,  and 
probability  modeling,  this  investigation  attempts  to  develop  a  model  to  predict  Tie  Bar  life. 
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n.  COMBINATORIAL  FAILURE  SEQUENCES 

O’Cotmor  [Ref.  1)  studied  the  Tie  Bar  problem  in  a  previous  investigation.  He  identified 
that  the  probability  of  failure  of  a  system  with  multiple  components  could  be  determined  in  terms 
of  the  probability  of  failure  of  the  components.  That  is,  the  probability  of  failure  of  the  Tie  Bar  can 
be  determined  in  terms  of  the  probability  of  failure  of  the  straps  which  can  in  turn  be  determined  in 
terms  of  the  probability  of  failure  sequences  of  the  components  within  the  straps.  He  also  identified 
that  each  component  of  a  Tie  Bar  strap  can  be  thought  of  as  a  chain  in  which  failure  of  the  weakest 
link  results  in  failure  of  that  entire  component.  This  lead  to  adoption  of  a  Weibull  distribution  for 
the  individual  strap  components. 


F{x)  =  1  -  exps 


Equation  II>1 


Applied  to  the  Tie  Bar, 


Equation  II-2 

where  (^)]  is  the  probability  of  failure  of  component  A  within  a  strap  under  stress  Xa  in 

the  i'*'  configuration  (one,  two,  three,  or  four  component  failures),  xa(L)  is  the  applied  stress  on 
component  A  caused  by  the  global  load  L,  p  is  a  scale  parameter  of  the  material  strength,  and  a  is 
a  shape  parameter  of  the  material  variability.  It  is  important  to  note  that  the  components  have 
identical  fiiilure  probability  parameters  (a,  P)  because  the  material  is  identical,  but  because  the 
stress  states  are  not  the  same,  each  component  within  a  strap  has  a  different  probability  of  failure. 
Finite  element  analysis  is  used  here  to  calculate  the  stress  states  for  input  into  the  probability  of 
failure.  Equation  II-2.  Because  of  the  local  stress  dependence  on  the  combinatorial  failure 
sequences,  the  probability  of  failure  of  the  Tie  Bar  is  no  longer  Weibull. 
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A  total  of  16  failure  states  exist  for  a  four  component  Tie  Bar  strap.  Failure  states  are 
numbered  using  a  binary  code.  A  1  represents  an  unbroken  Tie  Bar  strap  component  and  a  0 
represents  a  broken  component.  Table  II-l  shows  the  possible  failure  states  and  a  decimal 
equivalent  name  for  each  failure  state. 


Component  State 

A  B  C  D 

Decimal 

Equivalent 

1 

1 

1 

1 

_ 15_ 

1 

1 

i 

0 

14 

1 

1 

0 

1 

13 

1 

1 

0 

0 

12 

1 

0 

1 

1 

11 

1 

0 

1 

0 

10 

1 

0 

0 

1 

9 

1 

0 

0 

0 

8 

0 

1 

i 

1 

7 

0 

1 

1 

0 

6 

0 

1 

0 

1 

5 

0 

1 

0 

0 

4 

0 

0 

1 

1 

3 

0 

0  n 

1 

0 

2 

0 

0 

0 

1  i 

1 

0 

0 

0 

0 

0 

1  represents  an  unbroken  Tie  ] 

Bar  strap  component 

0  represents  a  broken  component 


Table  II-l  Tie  Bar  Strap  Failure  States 

O’Connor  [Ref  1]  developed  an  algorithm  to  calculate  the  probability  of  failure  of  a  Tie 
Bar  leaf  given  the  initial  configuration;  zero,  one,  two,  three,  or  four  components  failed  and  the 
final  configuration  one,  two,  three,  four,  or  ail  failed.  He  first  developed  an  algorithm  to  calculate  a 
three  component  system  and  found  that  the  local  stress  can  be  calculated  for  each  break 
configuration,  but  the  relative  strengths  cannot  be  determined  imless  the  feilure  path  which  lead  to 
that  configuration  is  known.  His  three  component  “failure  tree”  is  shown  in  Figure  II-l. 
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A=B=C  Equally  Weak 


B=C  Weakest 

- .  w 

_ M _ 

A=C  Weakest 


A=B  Weakest 


_A 

C 


Indicates 

Failure 


B  Weakest 

- » 


Figure  II-l  Three  Component  Failure  Tree 


The  algorithm  to  generate  the  probability  of  system  failure  given  the  component  failure 
first  requires  calculating  the  probability  of  failure  of  a  sequence  of  consecutive  states,  (one  branch 
or  path  of  the  failure  tree).  The  probability  of  failure  of  one  sequence  is  the  intersection  or  product 
of  the  transitional  probabilities.  The  probability  of  failure  of  the  system  is  the  union  or  sum  of  the 
probability  of  failure  of  each  of  the  sequences.  The  appendix  of  O’Connor’s  work  [Ref  1]  lists  all 
the  possible  failure  permutations  for  a  four  component  system  given  that  none  or  any  one,  two,  or 
three  of  the  four  components  has  already  failed. 
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in.  FINITE  ELEMENT  MODELS 

In  order  to  calculate  the  probability  of  failure  of  components  of  the  strap,  the  stress  state  at 
the  prospective  failure  location  must  be  known  for  each  of  the  physical  states,  (none  failed  and 
partial  failures  in  Figure  II-l).  Finite  element  analysis  is  used  to  calculate  such  states  for  three 
different  designs.  These  designs  include: 

A.  The  current  design 

B.  An  elliptical  slot  design  proposed  by  Boeing  Vertol 

C.  A  reduced  constraint  design  proposed  in  this  investigation 

A.  CURRENT  DESIGN 

A  Finite  Element  model  of  the  CH-46  Tie  Bar,  Figure  IIl-l,  was  created  using 
MSC/NASTRAN  for  Windows.  The  model  consists  of  the  top  and  bottom  straps  of  the  Tie  Bar, 
the  three  expanded  pins  at  each  end  of  the  Tie  Bar  that  hold  the  assembly  together,  and  the  two 
large  pitch  housing/pitch  shaft  pins  (parts  162  and  167  of  Figure  1-4)  that  connect  the  Tie  Bar  to 
the  pitch  shaft  and  pitch  housing.  The  straps  are  modeled  using  quad  PLATE  elements.  The  top 
and  bottom  straps  are  meshed  identically. 
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Figure  III-l  Current  Tie  Bar  Finite  Element  Model 


The  expanded  assembly  pins  and  pitch  housing/pitch  shaft  pins  are  modeled  using  RIGID 
elements.  The  assembly  pin  RIGID  elements  tie  all  degrees  of  freedom  of  a  bottom  strap  node  to  a 
corresponding  top  strap  node.  The  pitch  housing/pitch  shaft  pins  (Figure  III*2)  are  modeled  by 
two  vertical  RIGID  elements  to  which  “bicycle  spoke”  RIGID  elements  are  radially  attached. 
Boundary  conditions  are  applied  to  the  center  node  of  the  pitch  housing/pitch  shaft  pins  to  model 
the  displacements  of  the  Tie  Bar  in  flight. 


The  pitch  housing/pitch  shaft  pins  in  this  finite  element  model  are  attached  to  the  tie  bar 
straps.  These  pins  in  the  actual  Tie  Bar  are  not  attached.  Attaching  the  pins  to  the  Tie  Bar  straps 
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provides  expedient  modeling  of  the  equivalent  boundary  conditions.  The  local  stress  around  the  pin 
holes  is  affected  by  this  expediency;  however,  service  history  has  shown  that  failures  occur  in  the  4 
components  of  each  Tie  Bar  strap,  not  the  areas  around  the  holes.  Saint-Venant’s  principle  assures 
that  the  stresses  at  actual  failure  locations,  which  are  far  away  from  the  pin  holes,  are  accurately 
modeled.  By  applying  the  boundary  conditions  to  the  attached  model  pins,  all  the  displacement  is 
applied  directly  to  the  4  strap  components.  The  pitch  housing/pitch  shaft  pin  holes  remain  circular; 
they  do  not  elongate.  Therefore,  the  four  strap  components  experience  slightly  more  strain  than  the 
real  Tie  Bar,  and  the  model  is  conservative. 

The  Tie  Bar  straps  are  positioned  to  simulate  the  maximum  allowable  Tie  Bar  stack 
thickness.  The  PLATE  elements  around  the  pitch  housing/pitch  shaft  pins  are  thicker  than  the  rest 
of  the  Tie  Bar.  It  is  assumed  that  the  spacers  on  both  sides  of  the  outermost  Tie  Bar  straps 
contribute  to  the  strength  and  stiffiiess  of  that  region.  Figure  III-3  illustrates  the  thickened  region. 
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B.  ELLIPTICAL  SLOT  DESIGN 


The  Tie  Bar  manufacturer,  Boeing  Vertol,  has  suggested  modifying  the  geometry  of  the 
Tie  Bar  strap  slots  to  reduce  stress  concentration  and  increase  Tie  Bar  life.  The  focus  of  the 
modification  is  to  change  the  previously  semi-circular  slot  end  to  a  semi-elliptical  slot  end  as  shown 
in  Figure  111-4.  This  geometric  modification  was  made  to  the  finite  element  model  to  replicate  the 
new  geometry. 
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Current  and  Reduced  Constraint  Elliptical  Slot  Geometry 

Slot  Geometry 

Figure  111-4  Tie  Bar  Slot  Geometry  Modification 

C.  REDUCED  CONSTRAINT  DESIGN 

In  studying  the  kinematics  of  the  Tie  Bar,  it  was  observed  that  rotor  blade  feathering 
causes  straps  away  from  the  Tie  Bar  neutral  axis  to  experience  not  only  a  twisting  deformation  but 
also  a  parallel  shear.  If  the  three  swedged  assembly  pins  around  the  pitch  shaft/pitch  housing  pins 
are  removed,  the  reduced  constraint  significantly  minimizes  the  shear  deformation.  Without 
changing  the  geometry  from  the  current  Tie  Bar  design,  another  set  of  models  was  generated  by 
removing  the  three  swedged  assembly  pins  from  each  end  of  the  Tie  Bar  model. 
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Figure  I1I-5  Reduced  Constraint  Tie  Bar  Finite  Element  Model 


The  straps  for  the  reduced  constraint  models  are  only  held  together  by  the  pitch 
housing/pitch  shaft  pin  on  either  end.  In  a  real  Tie  Bar  of  this  design,  the  straps  could  be  held 
together  by  a  sleeve  with  rolled-over  ends.  Figure  111-6.  The  assembly  pin  holes  would  not  be 
required,  and  the  Tie  Bar  strap  ends  would  be  rounded  off  to  maintain  clearance  inside  the  pitch 


Figure  III-6  Sleeve  to  Assemble  Stack  of  Reduced  Constraint  Tie  Bar 

shaft  when  the  straps  rotate.  Also,  a  low  friction  material  such  as  Teflon  could  be  used  for  the 
spacers  between  straps  to  reduce  friction.  These  models  were  easy  to  create  from  the  current  Tie 
Bar  models  and  permitted  direct  comparison. 

Whereas  the  current  and  elliptical  slot  design  models  have  all  degrees  of  freedom  (DOF)  of 
the  pitch  housing/pitch  shaft  pin  end  nodes  tied  to  the  pitch  housing/pitch  shaft  pin  master  node, 
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the  reduced  constraint  model  allows  the  pitch  housing/pitch  shaft  pin  end  nodes  to  rotate  about  the 
pin  z-axis  relative  to  the  pitch  housing/pitch  shaft  pin  master  node.  This  design.  Figure  111-7, 
alleviates  in-plane  shear  and  bending  of  the  Tie  Bar  straps. 


Figure  III-7  Pitch  housing/pitch  shaft  Pin  with  Rotating  Ends 
D.  FAILURE  STATES 

Models  with  break  sites  were  created  from  the  unbroken,  “parent”  models.  Individual  top 
strap  components  were  “broken”  to  simulate  breakage  in  service.  This  was  accomplished  by 
disconnecting  the  nodes  at  the  prospective  ftiilure  sites  to  model  the  crack  trajectory.  With  the 
nodes  disconnected,  the  affected  PLATE  elements  in  the  model  were  shrunk  to  show  the  “break”  in 
the  Tie  Bar  strap. 

As  stated  previously,  a  total  of  16  ftiilure  states  exists  for  a  four  component  Tie  Bar  strap. 
The  possible  failure  states  are  shown  in  Table  II-l.  Six  of  these  states  are  eliminated  using 
symmetry,  and  the  all-broken  state  does  not  require  analysis.  So,  eight  models  were  “grown”  from 
the  original  model  to  recreate  the  possible  ftiilure  states.  Figure  III-8  illustrates  a  broken  Tie  Bar 
model. 
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Figure  III-8  Broken  Tie  Bar  Strap 


E.  BOUNDARY  CONDITIONS 

The  boundary  conditions  for  the  Tie  Bar  were  supplied  by  the  Tie  Bar  manufacturer  as  a 
flight  spectrum.  For  each  flight  condition,  values  of  centrifugal  force,  twist,  and  bending  were 
provided.  The  flight  conditions  consist  of  G-A-G,  Pre-Flight,  HSLF,  Maneuver,  and  Autorotation, 
and  they  are  listed  with  their  corresponding  Tie  Bar  displacements  in  Table  III- 1. 


Condition 


G-A-G 


Cycles 


Tie  Bar  Twist* 


0->113%->0  0->76250  lb->0  •25°->-16'’^+22°  0-)'±0.365“^0  1500 

0->298  RPM-»0 _ 

G-A-G  1  0^118%->0  0->83150  lb->0  -25°->-16°->+22°  0->±0.365°-+0  1500 


Pre-Flight 


Maintenance 


HSLF 


Maneuver  1 


0^312RPM->0 


0 


100%  Nr 
264  RPM 


105%  Nr 
277  RPM 


105%  Nr 
277  RPM 


105%  Nr 
277  RPM 


105%  Nr 
277  RPM 


113%  Nr 
298  RPM 


113%  Nr 
298  RPM 


113%  Nr 
298  RPM 


113%  Nr 
298  RPM 


11 8%  Nr 
312  RPM 


133%  Nr 
351  RPM 


59715  lb 


65835  lb 


65835  lb 


65835  lb 


65835  lb 


76250  lb 


76250  lb 


76250  lb 


76250  lb 


831501b 


105630  lb 


139960  lb 


-25°->+22' 

+50° 


0±6° 


Negligible 

0 


Once  per 
Lifetime 


N/A 


*  Twist  cycle  is  once  per  rotor  revolution  for  HSLF,  Maneuvers,  and  Autorotation  conditions. 

♦  Relative  angle  of  attachment  pins  referenced  to  the  vertical  pin.  Bending  cycle  is  two  times  per  rotor 
revolution. 


Table  Ill-l  Tie  Bar  Loading  Conditions 


1.  Flight  Condition  Description 

The  G-A-G  condition  is  a  ground-air-ground  condition.  The  first  half  of  the  G-A-G  cycle 
includes  engaging  the  aircraft  rotors,  checking  the  rotor  RPM  at  maximum  throttle  with  the 
collective  bottomed,  and  lifting  the  aircraft  off  the  ground.  The  second  half  of  the  G-A-G  cycle 
includes  the  aircraft  landing  and  rotor  shutdown. 
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Preflight  checks  occur  on  the  ground  with  the  rotors  disengaged.  The  collective  is  moved 
from  full  down  to  full  up.  The  cyclic  is  moved  full  forward,  full  aft,  full  left,  and  full  right  during 
these  checks.  These  control  movements  twist  the  Tie  Bar. 

HSLF  is  horizontal  steady  level  flight.  It  is  a  cruise  condition  without  climbs,  descents,  or 

turns. 

Autorotation  is  a  power  off  maneuver  in  which  the  collective  is  lowered  and  the  aircraft's 
descent  is  slowed  by  the  windmilling  of  the  rotor  blades.  As  shown  in  Table  III-2,  autorotation  is 
generally  flown  at  a  higher  rotor  RPM  than  the  other  flight  conditions. 

2.  Rotor  System  Behavior 

O’Connor  [Ref.  1]  conducted  a  kinematic  analysis  of  a  single  rotor  system  in  previous  Tie 
Bar  research  to  establish  the  phasing  of  the  loads  and  the  kinematics  of  the  Tie  Bar.  A  fully 
articulated  system,  the  CH-46  rotor  system  has  both  lead-lag  and  flapping  hinges  to  provide 
independent  in-plane  and  out  of  plane  motion  of  each  rotor  blade.  The  rotor  blade  flapping  hinge  is 
designed  into  the  rotor  system  to  prevent  development  of  an  undesired  rolling  moment  on  the 
helicopter.  Flapping  occurs  because  the  advancing  rotor  blade  sees  a  higher  relative  velocity  and 
hence  higher  lift  than  the  retreating  blade.  Maximum  upward  flapping  occurs  at  \|/  =  180°,  and 
minimum  flapping  occurs  at  vp  =  0°.  Through  his  research,  O'Connor  also  discovered  that  the 
maximum  flatwise  bending  moment  occurs  at  vp  =  90°— >120°  (upward)  and  vp  =  270°  (downward). 
The  rotor  angle  convention  is  shown  in  Figure  I1I-9. 


19 


VJ/  =  90®  V  =  270® 


Figure  III-9  Rotor  Angle  Convention 

Lead-lag  is  caused  by  cyclically  varying  drag  and  inertia  loads.  As  a  blade  flaps  up,  its 
center  of  gravity  moves  inward.  Since  angular  momentum  must  be  conserved,  the  blade  accelerates 
in  plane  to  compensate  for  the  decrease  in  moment  of  inertia.  The  change  in  inertia  of  this  upward 
flapping  blade  is  resisted  by  the  other  blades  and  an  in-plane  moment  is  produced.  The  lead-lag 
hinge  allows  the  blade  to  accelerate  forward.  As  the  blade  flaps  down,  the  center  of  gravity  moves 
out,  and  the  blade  decelerates. 

Maximum  twist  or  feathering  of  the  rotor  blade  occurs  at  vj/  =  270®  where  the  rotor  blade 
increases  pitch  to  compensate  for  the  lower  relative  velocity.  Minimum  feathering  occurs  at  vy  = 
90®  where  relative  velocity  is  highest. 
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O'Connor  [Ref.  1]  also  researched  the  phasing  of  the  Tie  Bar  bending  and  twist.  He  noted 
that  the  lead-lag  produced  negligible  bending  of  the  Tie  Bar.  His  study  of  the  rotor  blade  dynamics 
produced  the  results  in  Table  III-2. 

Rotor  Position  -  v(/  Centrifugal  Force  Feathering  Bending 

(degrees) 


0 

constant 

mean 

mean 

90 

constant 

minimum 

max  up 

180 

constant 

mean 

mean 

270 

constant 

maximum 

max  down 

Table  III- 2  Rotor  Phasing 


3.  Conversion  of  Flight  Conditions  to  Model  Boundary  Conditions 

Three  primary  load  cases  were  run  for  each  Tie  Bar  design  and  break  state.  The  three  load 
cases  were: 

1.  0.2  inch  extension  at  each  pitch  housing/pitch  shaft  pin.  (X-translation);  Figure  III-IO 

2.  0.2  radian  twist  at  each  pitch  housing/pitch  shaft  pin.  (X-rotation);  Figure  III-ll 

3.  0.2  radian  bending  at  each  pitch  housing/pitch  shaft  pin.  (Y -rotation);  Figure  III-12 
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Figure  III-IO  Primary  Extension  Load  Case  (Side  View) 


Figure  111-11  Primary  Twist  Load  Case 
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Figure  III-12  Primary  Bending  Load  Case 

These  three  basic  load  cases  were  scaled  and  linearly  combined  to  recreate  each  flight 
condition  at  vy  =  0°,  90°,  180°,  and  270°.  The  twist  and  bending  boundary  conditions  were  scaled 
directly.  The  extension  boundary  conditions  were  given  as  centrifugal  forces  instead  of 
displacements.  To  scale  these  boundary  conditions,  the  force  to  extend  each  unbroken  model  0.2 
inches  was  read  from  the  NASTRAN  output.  This  value  was  divided  by  2  to  find  the  force 
required  to  deform  1  unbroken  strap  0.2  inches.  This  result  was  multiplied  by  35  to  determine  the 
force  required  to  deform  a  whole  35-strap,  unbroken  Tie  Bar  0.2  inches. 

Equation  III-l 

where  force  required  to  deform  a  whole  Tie  Bar  0.2  inches,  and  (/^02)2,  force 

required  to  deform  a  two-strap  NASTRAN  model. 

To  calculate  the  extension  boundary  condition  scaling  fector  the  following  equation  was 

used: 
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CF 


Extension  _  scale  _  factor 


((F.2k-(F„)r,)  +  34*(F„);/ 


Equation  III-2 


where  CF  is  the  centrifugal  force  for  the  given  flight  condition,  and(FQ  2),^  is  the  force  required  to 
deform  one  unbroken  strap  0.2  inches. 
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IV.  FINITE  ELEMENT  MODEL  RESULTS 

Service  failures  of  the  current  Tie-Bar  and  finite  element  model  stress  contour  plots 
indicate  that  the  areas  around  the  Tie-Bar  slot  ends  are  most  critical.  At  every  spatial  location,  the 
3  scalar  components  of  the  stress  tensor  can  be  combined  to  form  a  scalar  value  based  on  the 
failure  process.  Since  stainless  steel  is  a  ductile  material,  the  von  Mises  failure  stress  is  used.  The 
NASTRAN  contour  plots  that  follow.  Figure  IV- 1,  Figure  IV-2,  and  Figure  IV-3  show  the 
deformation  and  von  Mises  stress  contour  of  the  top  strap  of  the  three  Tie-Bar  designs  when  the 
Maneuver  4  flight  condition  is  applied  at  vj/  =  270.  While  the  Maneuver  4  flight  condition  at  vj/  = 
270  loads  the  top  Tie  Bar  strap  more  heavily  than  the  bottom  one,  it  qualitatively  represents  typical 
behavior  for  the  models.  The  Current  Tie  Bar  and  the  Elliptical  Slot  Tie  Bar  models  exhibit  in¬ 
plane  bending  caused  by  the  assembly  pins.  The  assembly  pins  constrain  the  rotation  of  the  Tie  Bar 
straps  about  the  pitch  shaft/pitch  housing  pins.  The  Reduced  Constraint  Tie  Bar  model  exhibits  no 
in-plane  bending  since  the  assembly  pins  have  been  removed. 

The  current  Tie-Bar  model  shows  areas  of  high  stress  where  the  slots  transition  from  a 
straight  to  a  radius.  The  Elliptical  Slot  Tie-Bar  model  shows  stress  concentration  in  the  same 
areas,  but  at  a  slightly  lower  value.  The  Reduced  Constraint  Tie-Bar  Model  does  not  show  any  of 
the  stress  concentration  of  the  first  two  models  around  the  slot  ends. 

To  begin  understanding  the  effect  of  the  new  designs  on  the  stress  in  the  Tie-Bar 
quantitatively,  several  elements  of  interest  were  chosen.  These  elements  are  located  at  the  transition 
between  the  straight  part  of  the  slots  and  the  radiused  ends.  The  element  numbers  are  shown  in  the 
“element  map”.  Figure  IV-4. 
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Figure  IV-2  Elliptical  Slot  Tie-Bar,  Top  Strap,  Maneuver  4,  Y=270 


Figure  IV-3  Reduced  Constraint  Tie-Bar,  Top  Strap,  Maneuver  4,  Y=270 
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Figure  IV-4  Map  of  NASTRAN  Element  Numbers 


Focusing  on  these  “critical  elements”,  the  effect  of  the  design  changes  was  first  evaluated 
by  looking  at  the  “critical  elements”  with  the  primary  load  cases  applied  individually.  This  allowed 
study  of  individual  effects  on  the  Tie  Bar  behavior. 

For  the  extension  only  load  case,  the  Tie-Bar  was  elongated  0.02  inches,  a  typical  flight 
value.  When  loaded  in  extension  only,  the  stresses  for  the  Current  Tie-Bar  and  the  Reduced 
Constraint  Tie  Bar  are  identical  as  expected,  with  the  highest  values  occurring  at  the  inner  edges  of 
the  outer  slots.  The  maximum  stresses  for  the  Elliptical  Slot  model  seem  to  shift  away  from  the  slot 
end  slightly.  The  magnitude  is  not  always  lower.  This  suggests  that  the  elliptical  slot  design  may  be 
a  benefit  to  some  locations  and  a  detriment  to  others.  Figure  IV-5  summarizes  these  results. 


Extension  Only  (0.02  inches) 


a  Current  Tie  Bar  Model 
S  Elliptical  Slot  Tie  Bar  Model 
a  Reduced  Constraint  Tie  Bar  Model 


Nattnn  Element  Number 


Figure  IV-5  Comparison  of  Critical  Element  Stresses  for  Three  Designs,  Extension 

In  the  twist  only  load  case,  more  dramatic  results  are  evident.  Again,  the  Elliptical  Slot 
model  shifts  the  stress  slightly  away  from  the  slot  end.  It  also  decreases  the  maximum  von  Mises 
stress  by  nearly  10,000  psi  or  17  percent.  The  Reduced  Constraint  model  shows  even  more 
spectacular  results.  This  model  reduces  the  maximum  von  Mises  stress  by  3 1,000  psi  or  53 
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percent.  Not  only  does  the  Reduced  Constraint  model  drastically  reduce  the  maximum  von  Mises 
stress,  but  it  also  evens  out  the  stress  between  the  elements.  Figure  IV-6  summarizes  these  results. 


Twist  Only  (0.4  radians) 


a  Cunrent  Tie  Bar  Model 
B  Elliptical  Slot  Tie  Bar  Model 
B  Reduced  Constraint  Tie  Bar  Model 


267  268  271  272  287  288  289  290  303  304  305  306  307  308  317  318 

Nastnn  Element  Number 


Figure  IV-6  Comparison  of  Critical  Element  Stresses  for  Three  Designs,  Twist 

The  bending  only  load  case  causes  the  models  to  exhibit  behavior  similar  to  that  produced 
by  the  extension  load  case.  Again,  the  stress  for  the  Elliptical  Slot  model  moves  slightly  away  from 
the  slot  ends.  The  stresses  produced  by  the  bending  load  case  are  only  about  one-fifth  those 
produced  by  extension  and  twist.  Figure  IV-7  summarizes  these  results. 
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Figure  IV-7  Comparison  of  Critical  Element  Stresses  for  Three  Designs,  Bending 

To  summarize,  the  stresses  due  to  bending  are  only  one-fifth  of  those  due  to  extension  and 
twisting.  The  bending  mode  is  judged  to  be  not  critical.  Extension  induced  stresses  at  all  critical 
locations  are  within  a  narrow  range.  The  elliptical  slot  mainly  shifts  the  maximum  stress  locations; 
its  benefit  is  not  definitely  obvious.  The  twisting  induced  stresses  are  dramatically  lowered  by  50 
percent  by  the  reduced  constraint  configuration.  Serious  consideration  of  this  modification  is 
warranted. 

Results  of  the  Maneuver  4  flight  load  case  at  v|/  =  270  are  shown  in  Figure  IV-8.  These 
results  are  formed  from  a  scaled  linear  combination  of  the  above  uncoupled  load  cases.  This 
combmed  load  case  indicates  that  the  Elliptical  Slot  design  reduces  the  maximum  stress  by  about 
5,000  psi  or  5  percent.  The  Reduced  Constraint  design  reduces  maximum  stress  by  about  25,000 
psi  or  26  percent.  Again  it  is  apparent  that  the  Reduced  Constraint  design  evens  out  the  stresses 
among  the  elements.  In  feet,  the  Reduced  Constraint  design  shows  an  increase  in  stress  for  some 
elements,  but  the  maximum  stress  decreases.  It  is  important  to  remember  that  these  results  are  for 
only  one  flight  condition  at  one  rotor  position.  The  trend,  however,  is  significant. 
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Maneuver  4,  ‘F  =  270* 


a  Current  Tie  Bar  Model 
a  Elliptical  Slot  Tie  Bar  Model 
B  Reduced  Constraint  Tie  Bar  Model 


100.000 
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80,000 

70,000 

60,000 
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30,000 

20,000 
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267  268  271  272  287  288  28®  290  303  304  306  308  307  308  317  318 

Nactrtn  Element  Number 


Figure  IV-8  Comparison  of  Critical  Element  Stresses  for  Three  Designs,  Maneuver  4 

In  order  to  evaluate  the  cumulative  effect  of  the  three  designs  at  the  different  rotor 
orientations,  results  were  extracted  from  the  finite  element  analysis.  Figure  IV-9  and  Figure  IV-10 
show  how  the  von  Mises  stress  varies  for  one  chosen  element  as  the  rotor  moves  through  one 
revolution.  The  element  chosen  was  the  maximum  stress  element  for  the  flight  condition.  The  solid 
lines  are  the  current  design  stresses,  the  dotted  lines  are  the  elliptical  slot  design  stresses,  and  the 
dashed  lines  are  the  reduced  constraint  design  stresses.  The  elliptical  slot  design  lowers  both  the 
mean  stress  and  stress  amplitude  slightly.  The  reduced  constraint  design  lowers  both  the  mean 
stress  and  stress  amplitude  significantly  for  the  Maneuver  4  condition  and  lowers  the  mean  stress 
for  the  Autorotation  5  condition. 
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von  Misos  Stress  (psi)  , — ,  Mises  Stress  (psi) 


Maneuver  4 


Autorotation  5 


Figure  IV-10  Slot  End  Individual  Element  von  Mises  Stresses,  Autorotation  5 
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It  is  also  interesting  to  note  that  the  moment  required  to  twist  the  reduced  constraint  design 
is  significantly  lower  than  that  required  for  both  the  current  and  elliptical  slot  designs.  The  moment 
required  to  twist  the  modified  slot  design  is  greatest.  The  NASTRAN  R1  Constraint  Moments  for 
the  three  designs  with  a  0.4  radian  twist  applied  are  tabulated  in  Table  IV-1. 


Tie  Bar  Design 

R1  (Twist)  Constraint  Moment* 

(in-lb) 

Current 

43.3 

Elliptical  Slot 

44.0 

Reduced  Constraint 

13.0 

♦0.4  radian  twist  of  the  two-strap  NASTRAN  model 

Table  IV-1  Tie  Bar  Constraint  Moments 
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V.  LIFE  PREDICTION 


To  effectively  predict  Tie  Bar  life  and  evaluate  design  changes,  the  probability  density 
function,  PDF,  and  the  cumulative  distribution  function,  CDF,  for  the  tie  bar  system  are  needed. 
Determination  of  the  forms  of  these  distributions  requires  partitioning  the  Tie  Bar  system  mto  its 
components  whose  distribution  functions  are  known.  Conceptually,  each  Tie  Bar  strap  is  first 
partitioned  into  4  components.  Figure  1-6.  Each  of  these  components  is  in  turn  partitioned  into 
smaller  volumes  or  finite  elements.  In  a  finite-element  model  with  all  elements  of  comparable  size, 
any  element  could  describe  a  metric  volume.  The  size  of  this  metric  is  arbitrary;  a  larger  metric  can 
be  thought  of  as  a  combination  of  smaller  metric  volumes.  A  larger  metric,  however,  has  a  higher 
probability  of  failure  and  a  lower  strength.  The  selection  of  the  metric  volume  is  such  that  the 
random  variable  is  spatially  homogeneous;  that  is,  the  stresses  are  uniformly  distributed.  The 
uniform  stress  criterion  is  the  same  as  selection  for  element  size  in  finite  element  analysis;  as  a 
result,  the  element  can  be  used  directly  as  a  metric. 

A.  FLAW  DISTRIBUTION 

The  mechanistic  bases  for  failure  in  metals  are  by  flaw  growth  and  by  dislocations.  For 
brittle  failure  modes,  a  given  element  contains  flaws  that  are  either  smaller  than  some  critical  flaw 
size,  a  <  a*,  or  larger,  a  >  a*.  Each  element  can  be  considered  flawed  or  not  flawed  based  on  this 
criterion.  In  other  words,  scanning  each  element  for  flaw  size  would  create  two  categories  of 
elements,  those  with  flaw  sizes  greater  than  or  equal  to  critical,  and  those  with  flaw  sizes  less  than 
critical.  There  are  only  two  possibilities  for  each  element,  pass  or  fail.  Respective  to  the  critical 
flaw  size,  the  exceedences  are  binomially  distributed.  The  binomial  distribution  has  the  form: 

Equation  V-1 

where  f(n)  is  the  probability  density  function,  pdf. 
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"  iN-n)\n\  ’ 

Equation  V-2 

p  is  the  probability  that  the  element  has  a  >  ac,  n  is  the  number  of  elements  that  have  a  >  a^,  and 

N  is  the  total  number  of  metric  elements  being  considered.  In  the  current  investigation,  it  is  the  total 
number  of  elements  in  each  strap. 

When  the  number  of  elements  is  large,  and  the  flaw  density  is  low,  then  the  probability 
density  function,  pdf,  takes  on  the  Poisson  distribution  form: 


Equation  V-3 

where  psNp  and  is  assumed  to  remain  constant. 

For  fatigue,  the  location  parameter  of  the  distribution  of  flaws  is  time  dependent.  As  time 
increases,  flaws  grow  and  the  probability  that  flaws  exceed  the  critical  size  increases.  The  binomial 
distribution  then,  is  a  snapshot  in  time. 

B.  COMPONENT  LIFE  DISTRIBUTION 

For  any  given  instant  in  time,  in  order  for  a  given  component  to  have  life  x,  each  of  its 
elements  must  have  life  x.  In  other  words,  element  I  and  element  2  and  element  3  and  ...  element  n 
must  have  life  x  in  order  that  the  entire  system  has  life  x.  Conversely,  if  any  one  element  does  not 
have  life  x,  the  whole  component  foils.  In  other  words,  when  the  weakest-link  foils,  the  system  fails. 
In  terms  of  the  previous  discussion,  when  any  element  experiences  crack  growth  to  a  size  greater 
than  the  critical  dimension,  the  whole  Tie  Bar  component  fails. 

In  Boolean  algebra,  the  weakest  link  system  reliability  is  the  intersection  (n).  The 
reliability  of  the  components  is  conceptually  a  chain,  or  series,  of  Imks  with  no  parallel  or 
redundant  links.  For  a  three  component  system  where  R  denotes  reliability  or  probability  of  non¬ 
failure. 
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R,{A-,nA-,nA-,}  = 

Equation  V-4 

Assuming  now  that  the  component  is  made  up  of  many  elements  and  the  reliability  of  each  is 
independent. 

Equation  V-5 

where  R„  denotes  the  Tie  Bar  component  reliability.  In  this  case,  the  component  being  discussed  is 
itself  a  system  of  elements  or  links  in  the  chain.  Failure  is  assumed  homogeneous  to  mechanism.  In 
other  words,  failure  is  not  related  to  the  size  of  a  component.  A  larger  component  will  have  more 
elements,  but  the  flaw  distribution  does  not  change.  When  this  is  true,  a  hazard  function  ^(t)  is 
defined  which  takes  on  the  Weibull  form. 


Eir)^r\ 

Equation  V-6 

where  H(t)  is  the  system  reliability  or  hazard  function,  x  is  the  life  of  a  single  element,  and  a  is  the 
number  of  elements  in  the  component.  The  cumulative  distribution  function,  CDF,  can  then  be 
written  as: 

F(0  =  l-exp[-5(r)]. 

Equation  V-7 

It  is  important  to  recognize  that  the  Weibull  life  distribution  contains  an  underlying  Poisson  flaw 
distribution,  which  was  described  in  the  previous  sub-section. 

C.  LIFE  CONVOLUTION 

The  intrinsic  normalized  life,  x,  for  a  given  stress  history,  S(t),  is  obtained  by  convoluting  the  effect 
of  stress  via  the  breakdown  rule  k. 
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Equation  V-8 

where  S(t)  is  a  stress  norm  that  is  piecewise  continuous  in  time  t,  f  is  a  non-dimensionalizing  and 
normalizing  parameter  for  time,  ti  is  the  initial  time,  tf  is  the  final  time,  and  k  is  a  damage  function. 
In  the  simplest  case,  this  process  reduces  to  the  familiar  Miner’s  Rule,  which  is  a  linear 
superposition  of  damage. 

Different  physical  processes  give  rise  to  different  forms  of  the  damage  fimction  K{S{t)) . 
The  power  form  and  exponential  form  are  explored.  Combinations  of  the  two  forms  are  also 
possible. 


1.  Power  Law  Damage  Function 

The  first  proposed  damage  function  is  defined  using  a  power  relation.  This  form  has  been 
observed  to  fit  high  cycle  fatigue  data  in  metals  and  is  associated  with  flaw  growth. 

Equation  V-9 

where  b  is  a  constant  exponent,  and  C,  is  a  constant  non-dimensionalizing  parameter  for  stress. 

Both  b  and  Cl  are  material  constants.  The  constants  are  determined  by  fitting  a  Ime  to  Strength- 
Life  data  for  a  given  material  in  log  space  where 

1 

slope 

Equation  V-10 
and 
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Cl  =  exp(intercept). 


Equation  V-11 

Substituting  Equation  V-9  into  Equation  V-8  yields, 


Equation  V-12 

and  substituting  Equation  V-12  into  Equation  V-6  yields. 


Equation  V-13 


Finally,  substituting  Equation  V-13  into  Equation  V-7  yields; 


F(t\S)  =  1  -  exp 


Equation  V-14 

Next,  the  parameters  from  the  life  test  are  related  to  the  parameters  in  standard  Weibull 
form.  Standard  Weibull  form  is: 


i?(0  =  l-F(0  =  exp 


“f  ' 

t 

Equation  V-15 

So,  Equation  V-14  is  rearranged  and  equated  to  Equation  V-15. 
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Equation  V-16 

Every  element  in  Equation  V-6  is  now  known.  /  is  an  arbitrary  normalizing  parameter. 
Cl,  b,  tp,  and  Ot  are  determined  from  material  testing,  and  a  is  the  size  effect  parameter,  which  is 
unity  if  the  test  specimens  are  the  same  size  as  the  actual  part.  Given  a  stress  history.  Equation  V- 
6  can  be  solved  for  life  t. 


2.  Exponential  Law  Damage  Function 

The  second  proposed  damage  function  is  defined  using  an  exponential  law.  This  form  has 
been  observed  to  fit  low  cycle  fatigue  data  in  metals  and  is  associated  with  yielding. 


K{S(t))  =  Trexp 
'^2 


Equation  V-17 

where  Cz  is  a  constant,  and  C3  is  a  constant  non-dimensionalizing  parameter  for  stress. 

Both  Cl  and  C3  are  material  constants.  The  constants  are  determined  by  fitting  a  line  to  Strength- 
Life  data  for  a  given  material  in  semi-log  space  where: 

C3  =  slope , 


and 


Equation  V-18 

Cl  =  exp(intercept  /  C3), 

Equation  V-19 


Substituting  Equation  V-17  into  Equation  V-8  yields: 
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Equation  V-20 

and  substituting  Equation  V-20  into  Equation  V-6  yields: 


Equation  V-21 

Substituting  Equation  V-21  into  Equation  V-7  yields: 


I  1  fv  1  (  S{t)  I  , 

F(,)=l-exp  ^expl— 


Equation  V-22 

Next,  the  parameters  from  the  life  test  are  related  to  the  parameters  in  standard  Weibull  form. 
So,  Equation  V-22  is  equated  to  Equation  V-IS. 


1  rv  1  (  5(0  L  ‘ 

exp  -  ;  J  Trexp  —  dt  =  exp  -  - 
\t  ’'<1  C,  \  L,  J  j 


Equation  V-23 


and  the  equation  is  solved  for  t. 


3.  Constant  Stress  Case  Example 

Having  developed  the  above,  the  solution  approach  for  the  power  law  form  of  the  damage 
function  is  to  determine  the  parameters  a,  b,  and  Ci  from  a  known  stress  history.  Once  known, 
these  parameters  can  be  used  to  predict  F(t)  of  another  history. 

For  example,  for  a  constant  stress  case  defined  as. 


S{t)  = 


0  t<t, 

s,  t>t, 


Equation  V-24 
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F{t)  =  1  -  exp 


CJ  V  / 


Equation  V-28 

which  is  the  cumulative  distribution  function  or  cumulative  probability  of  failure  for  the  constant 
stress  history. 

Next,  the  parameters  from  the  life  test  are  related  to  the  parameters  in  standard  Weibull 
form.  Equating  Equation  V-28  with  Equation  V-15, 


exp  —  — —  =exp  -  — 

^C,j  V  t  J  {ts 


Equation  V-29 


and  the  equation  is  solved  for  t. 


4.  Varying  Stress  Case  Based  on  Rotor  Kinematics 
The  Tie  Bar  stress  history  is  taken  to  be  of  the  form: 

SiHO)  =  S„-S^  cosiv^iO  -  90') , 


where 


Equation  V-30 


Equation  V-31 


-5. 


Equation  V-32 
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von  Mises  Stress  (psi) 


and  ^it)  is  the  rotor  position  in  degrees.  Figure  V-2  shows  a  typical  stress  history  for  the  Tie  Bar. 
Substituting  this  new  stress  history  into  Equation  V-12, 


Equation  V-33 

Unlike  the  previous  constant  load  case,  the  integrand  in  Equation  V-33  is  not  symbolically 
integrable.  Therefore,  the  integral  must  either  be  solved  numerically  or  approximated. 

Tie  Bar  Stress  History 


Rotor  Position,  'V  (degrees) 


Figure  V-2  Typical  Element  Stress  History 
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VI.  COUPON  TESTS 


Coupon  tests  were  performed  by  the  Tie  Bar  manufacturer  to  evaluate  different 
manufacturing  processes.  The  test  coupons  were  single,  modified  production  Tie-Bar  straps.  Each 
coupon  had  the  center  two  components  removed  by  machining. 


Figure  VI-1  Tie-Bar  Coupon 

The  following  manufacturing  process  modifications  were  evaluated  either  individually  or 
in  combination: 

•  Stamping  from  new  die 

•  Pre-twisting  to  1800  micro-strain 

•  Production  glass  bead  peening 

•  High  intensity  peening 

•  Peening  at  a  45  degree  angle 

•  Chemical  milling 

•  Deburring  and  steel  shot  peening 

A.  STAMPING  FROM  NEW  DIE 

A  new  die  was  used  to  stamp  some  of  the  coupons  to  get  a  cleaner  edge.  The  original  die 
creates  a  sheared  edge  that  is  fairly  rough  and  is  thought  to  provide  numerous  crack-initiation  sites. 

B.  PRE-TWISTING  TO  1800  MICRO-STRAIN 

Control  system  checkouts  are  performed  on  the  ground  with  the  rotors  stationary.  These 
control  checks  twist  the  Tie-Bar  approximately  22.5®.  Some  bending  also  occurs  due  to  the  rotor 
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blade  weight.  A  strain-gaged  Tie-Bar  installed  on  an  aircraft  showed  a  maximum  compression 
strain  of  -1800  micro-inches  during  this  control  check. 

To  emulate  this  check,  some  coupons  (two  at  a  time)  were  put  in  a  twist  test  fixture  with  a 
spacer  between  them  so  they  would  experience  the  same  loading  as  the  outermost  straps  of  an 
entire  Tie-Bar  assembly.  The  coupons  were  twisted  to  -1800  micro-inches  as  read  from  a  strain 
gage  mounted  on  the  outer  strap  near  a  slot  end  .  The  coupons  were  then  removed  from  the  twist 
fixture  and  tested  in  simple  tension-tension  until  failure. 

C.  6°  TWIST  LOAD 

Some  specimens  were  subjected  to  an  added  load  calculated  by  twisting  the  coupon  6  . 
Each  coupon  was  not  actually  twisted.  Instead,  a  strain  gage  was  mounted  on  a  coupon  at  the  slot 
radius  and  the  coupon  was  mounted  in  the  twist  test  jig.  The  coupon  was  twisted  6®  and  the  strain 
was  measured.  This  strain  was  converted  to  stress.  The  additional  coupon  load  was  calculated 
assuming  this  stress  equated  to  a  nominal  (P/A)  stress.  The  coupon  was  then  subjected  to  this  twist 
load  plus  the  standard  fatigue  load. 

The  coupons  were  tested  in  a  tension-tension  fetigue  machine  with  a  sinusoidal  loading. 
The  load  blocks  were  repeated  until  500  blocks  were  complete  or  the  coupon  failed,  except  as 
noted.  The  coupon  load  blocks  are  listed  in  the  appendix  Table  VIII- 1.  A  table  of  the  various 
configurations  and  their  associated  lives  are  also  located  in  the  appendix  Table  VIII-2. 

D.  DETERMINING  MATERIAL  PARAMETERS 

The  coupon  geometry  was  modeled  using  two  fimte-element  models:  one  for  the  standard 
coupons,  and  one  for  the  chemically  milled  coupons.  The  chemically  milled  coupons  had  0.002 
inches  of  material  removed  from  each  surface,  and  the  second  finite  element  model  reflects  this. 
The  basic  load  case  for  the  models  was  a  0. 1  inch  displacement  on  each  end  of  the  coupon  to  put 
the  coupon  in  tension. 
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Figure  VI-2  Coupon  Finite  Element  Model 


The  T1  Constraint  Force  was  calculated  for  each  model  and  was  scaled  for  each  coupon 
test  Load  Profile.  These  scale  factors  were  applied  to  the  fimte  element  models  to  recreate  the 
coupon  load  conditions.  Maximum  coupon  von  Mises  stresses  were  recorded  for  each  Load  Profile. 
These  maximum  von  Mises  stresses  were  converted  to  alternating  stresses,  S,  using  the  Goodman 
relation. 


Equation  VI-1 

where  Sn  is  the  fully  reversed  stress  level  corresponding  to  the  same  life  as  that  obtained  with  the 
stress  conditions,  a,  and  ffm,  and  o,  and  are  the  von  Mises  stress  amplitude  and  mean  stress 
determined  from  the  finite  element  model.  The  results  are  tabulated  in  the  appendix  and  are  plotted 
in  Figures  V-3  to  V-10.  These  figures  have  been  fitted  with  lines  using  a  least  squares  method.  In 
the  least  squares  method,  each  data  point  is  equally  weighted  and  is  an  expedient  method  of 
determining  the  central  tendency  of  the  fatigue  behavior,  i.e.  the  mean  of  the  S-N  curve. 
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In  Alternating  Stress,  von  Mises 


—•Linear  (Prcxluction  Part,  Produrtion  Peen) 


Figure  VI-3  Log-Log  S-n  Plot  of  Coupon  Test  Data,  Production  Part  and  Peen 


a  New  Stamping 
—Linear  (New  Stamping) 


Figure  VI-4  Log-Log  S-n  Plot  of  Coupon  Test  Data,  New  Stamping 
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8  9  10  11 

X  New  Stamping,  HMntensityPeen 
—Linear  (New  Stamping,  Hi-lntensity  Peen) 


Figure  VI-6  Log-Log  S-n  Plot  of  Coupon  T( 
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In  Alt«rnating  Stress,  von  Misss 


9  10  11  12  13  14  15  16 

•  New  Stamping,  Che»n-mill8d  ^ 

—Linear  (New  Stamping,  Chem-miHed) 


Figure  VI-7  Log-Log  S-n  Plot  of  Coupon  Test  Data,  New  Stamping,  Chem-Miiled 


+  New  Stamping.  Production  Peen,  Chem-miHed  ’fi  ^ 

—“Linear  (New Stamping,  Production  Peen,  Chem-miHed) 

Figure  VI-8  Log-Log  S-n  Plot  of  Coupon  Test  Data,  New  Stamping,  Production  Peen, 
Chem-Milted 
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Stress,  von  Mises 


•““Linear  (New  Stamping,  Oeburred) 


Figure  VI-9  Log-Log  S-n  Plot  of  Coupon  Test  Data,  New  Stamping,  Deburred 


“—Linear  (Naw  Stamping,  Chem-milled,  BlipticaO 


Figure  VI-10  Log-Log  S-n  Plot  of  Coupon  Test  Data,  New  Stamping,  Chem-Miiled, 
Elliptical  Slot 
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For  each  alternating  stress  level  in  Figures  V-3  to  V-10,  there  are  several  values  of  life. 
There  is  a  mean  distribution  for  each  of  these  stress  levels.  The  fitted  lines  can  be  interpreted  as  the 
mean  life  of  the  population  at  each  stress  level  and  may  be  corrected  to  the  location  parameter  by 
the  relation 

Equation  Vl-2 

For  this  conversion,  the  shape  parameter  a  must  be  known.  To  estimate  a,  a  maximum  likelihood 
estimator  technique  can  be  used  to  find  the  parameters  of  a  Weibull  distribution  that  fits  the  data  at 
each  stress  level.  Such  extensive  data  interpretation  has  not  been  pursued  in  the  investigation. 
Instead,  it  was  assumed  that  the  strength-life  relation  is  homologous;  that  is,  the  fractional  life  S-N 
curves  do  not  criss-cross.  Under  this  assumption,  the  shape  parameter  for  strength  oc,  can  be 
related  to  the  shape  parameter  for  life  cq  by 

b 

a,=  — , 

Equation  VI-3 

where  b  has  already  been  determined  by  the  least  squares  fit.  Typically,  strength  data  are  more 
plentiful  and  the  coefficient  of  variability  cv„  in  strength  are  reported  in  open  literature  to  be  0.5%. 
The  coefficient  of  variation  can  be  converted  to  the  shape  parameter  by  the  approximation 

1.2 

Cin  =  - • 

cv. 

Equation  Vl-4 

Table  VI-1  tabulates  the  material  parameters  for  the  power  law  development  in  the 
previous  section.  Larger  values  of  b  indicate  longer  life. 
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Condition 

b 

Ci  (psi) 

Production  Part,  Production  Peen 

-6.16 

364,171 

New  Stamping 

-2.77 

2,121,151 

Production  Part,  6°  Twist 

-7.26 

237,081 

New  Stamping,  Hi-intensity  Peen 

-13.61 

143,599 

New  Stamping,  Chem-milled 

-19.23 

168,054 

New  Stamping,  Production  Peen,  Chem-milled 

-13.21 

195,128 

New  Stamping,  Deburred 

-14.03 

123,731 

New  Stamping,  Chem-milled,  Elliptical 

-27.36 

172,509 

Table  Vl-1  Coupon  Test  Calculated  Parameters 
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vn.  CONCLUSIONS 


This  research  developed  a  probability  modeling  method  that  can  aid  the  prediction  of  a  Tie 
Bar  design  to  meet  desired  life  and  reliability  levels.  The  model  allows  the  effects  of  changes  in 
geometry,  material,  and  manufacturing  processes  to  be  quantified.  Since  the  model  can  predict  life, 
it  could  be  used  for  the  current  Tie  Bar  to  help  determine  a  safe  time  interval  between  inspections. 

Two  design  concepts  were  explored  in  this  investigation.  The  elliptical  slot  design  did  not 
reduce  stress  concentration  in  all  flight  kinematic  configurations.  Whether  the  elliptical  slot  design 
will  yield  longer  life,  is  not  definitive.  The  reduced  constraint  design  proposed  herein  provided 
much  more  stress  concentration  reduction  as  suggested  by  the  finite  element  study.  This  reduced 
constraint  design  may  produce  large  gains  in  life  over  the  current  design. 


57 


58 


vra.  RECOMMENDATIONS 


Some  field  data  has  been  recorded  on  current-design  Tie  Bars.  It  is  recommended  that 
detailed  breakage  data  including  part  hours,  number  of  broken  components,  and  location  of  broken 
components  continue  to  be  recorded.  This  data  can  be  used  to  calibrate  the  model  created  in  this 
investigation. 

While  the  modeling  method  developed  is  basically  complete  and  parts  of  it  have  been  used 
to  calculate  results,  more  effort  needs  to  be  spent  on  organizing  the  NASTRAN  output  and  feeding 
it  through  the  convolution  and  probability  models.  The  coupon  data  could  be  analyzed  more 
thoroughly  to  provide  better  material  parameters.  This  investigation  only  focused  on  one  strap  (4 
components)  of  a  two-strap  model.  Further  work  should  be  done  to  predict  the  probability  of 
failure  of  at  least  the  eight  components  of  the  two  outer  straps. 

Boeing  Vertol  is  qualifying  new  designs  for  the  Tie  Bar  through  testing.  Results  of  this 
testing  need  to  be  compared  to  the  model  results  to  ensure  the  model's  validity.  Field  data  can  be 
compared  to  both  the  model  results  and  the  testing  results  to  determine  what  further  load  and  life 
modeling  needs  to  be  done. 

Better  automation  and  generalization  of  the  methods  described  in  this  investigation  would 
be  quite  useful.  A  finite  element  post-processor  could  theoretically  be  developed  to  study  the  Tie 
Bar  as  well  as  any  general  structure. 
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APPENDIX 


*  Applied  One  Time  Every  Five  Blocks 
‘  1,000,000  cycles 
^  10,000,000  cycles 

Alternating  Stress  is  calculated  using  a  Croodman  relationship  and  von  Mises  stress  from 
the  coupon  finite  element  model. 

Table  VIII-1  Coupon  Load  Conditions 
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Specimen 

ID 

Production 

Part 

New 

Stamping 

6* 

Twist 

Production 

Peen 

High- 

Intensity 

Peen 

45’ 

Peen 

Chem- 

milled 

Debarred 

Elliptical 

Load 

Profile 

Life 

(cycles) 

45 

X 

X 

R 

39.250 

46 

X 

X 

C 

2.550.000 

47 

X 

X 

C 

48 

X 

X 

X 

2.270, 000| 

49 

X 

■■ 

X 

Table  VUI-2  Coupon  Test  Summary 
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